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^ ■ Abstract 



The underdamped Langevin equation of motion of a particle, in a symmetric periodic potential 



i-rt ! and subjected to a symmetric periodic forcing with mean zero over a period, with nonuniform 

Q ' friction, is solved numerically. The particle is shown to acquire a steady state mean velocity at 



asymptotically large time scales. This net particle velocity or the ratchet current is obtained in 



<^ , a range of forcing amplitudes Fq and peaks at some value of Fq within the range depending on 

Q"^ \ the value of the average friction coefficient and temperature of the medium. At these large time 

Tlj- ' scales the position dispersion grows proportionally with time, t, allowing for calculating the steady 

;<— ^ . state diffusion coefficient D which, interestingly, shows a peaking behaviour around the same Fq. 

OO ! 

(•~^ , The ratchet current, however, turns out to be largely coherent. At intermediate time scales, which 

bridge the small time scale behaviour of dispersion~ t^ to the large time one, the system shows, in 

j_j ■ some cases, periodic oscillation between dispersionless and steeply growing dispersion depending 

on the frequency of the forcing. The contribution of these different dispersion regimes to ratchet 

current is analysed. 
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I. INTRODUCTION 

The investigation of particle motion in periodic potentials has obvious relevance in con- 
densed matter studies. Motion of ions in a crystalline lattice is a case in point. Stochasticity 
in the motion is naturally introduced at nonzero temperatures. In these environment the par- 
ticle motion can be approximately described by a Langevin equation with suitable model po- 
tentials. Depending on the problem at hand the motion is either considered heavily damped, 
almost undamped, or in the intermediate situation mildly damped (or underdamped) . In 
many a situations in the former two extreme cases the Langevin equation becomes amenable 
to analytical solution. However, in the underdamped situation, barring a few special cases. 



numerical methods are used to solve the equation of motion of the particle 



l|. Owing to 



various kinds of errors and approximations involved in these (numerical) methods exact 
quantitative solutions are not possible. However, the method can reveal useful qualitative 
trends in the behaviour of the particle motion. For instance, recently it was shown'Z] that 
a particle, moving in an asymmetric but periodic potential in the presence of thermal noise 
when subjected to a symmetric periodic external drive (which adds to zero when averaged 
over a period), acquires a net motion when the parameters of the problem are chosen suit- 
ably. Such a net particle current without the application of any external bias or potential 
gradient in the presence of thermal noise is called thermal ratchet current and the system 
giving such a current is termed as thermal ratchet 3[] . Here the equilibrium condition of 
detailed balance is not applicable because the system was driven far away from equilibrium 
by rocking it periodically in the presence of noise. It was further shown that this system 



ute negative mobility 



J] . This prediction has already been found to be 



can even exhibit abso 

true experimentally|5|. It shows that in underdamped conditions or in the inertial regime 
diverse possibilities can be (qualitatively) uncovered by (numerically) solving the equations 
of motion. 

In the above important example the potential asymmetry was one of the necessary con- 
ditions for realization of ratchet current. The particle had to surmount the same potential 
barrier on either direction; only the slopes leading to the top of the barrier differed. A 
sinusoidal potential, for example, having no such asymmetry would not have yielded the 
ratchet current. In the present work, we consider similar particle motion in a sinusoidal 
potential. However, instead of a uniform friction coefficient of the medium we consider a 



model nonuniform space-dependent friction coefficient 7(x) of the medium. In particular, 
we consider a sinusoidally varying 7(2;) exactly similar to the potential but with a phase lag, 
0. A simple illustrative example of the model can be imagined thus: a stationary pressure 
wave is established in air giving a periodic 7(2;) for particle motion along x. An array of ions 
with the periodicity of 7(0;) but shifted a little to give a phase lag (f> will just fit our model 
for a charged particle motion along x. Here the potential is symmetric and periodic. How- 
ever, the directional symmetry of the system is disturbed by a phase shift in the similarly 
periodic 7(x). Particle motion along an one dimensional semiconductor heterostructure or 
protein motor motion on the surface of a microtubule along its axis would be some prac- 
tical situations close to the above example. Though we do not have any microscopic basis 
for justifying the model form of 7(2;) a periodic variation of friction has been argued ear- 
lier from mode-coupling theory of adatom motion on the surface of a crystal of identical 



atoms 



0. 



Also, the equation of motion has a direct correspondence with the resistively and 



capacitatively shunted junction (RCSJ) model of Josephson junctions; the term describing 
the nonuniformity of friction having an one-to-one correspondence with the 'cos (/>' term in 
the RCSJ model [7i]. A qualitative physical argument for the possibility of obtaining ratchet 
current in this inhomogeneous system was given in an earlier work[8]. 

We drive the system with a square-wave periodic field. The resulting Langevin equation 
is solved numerically. We obtain particle current and properties associated with it in the 
parameter space of external field amplitude Fq, the average friction coefficient 70, the phase 
lag (f), and the temperature T. Since it is a formidable task to explore the entire param- 
eter space, we present results for only some regions of a few sections of this space where 
appreciable ratchet current is obtained. 

The ratchet current is obtained in the steady state situation which is achieved in the 
asymptotic time limit. In our case we observe particle motion for a long time t such that the 
position dispersion ((Ax(t))^) averaged over many similar trajectories reach the situation 
where {{/\x{t))'^) ~ t. If the phase lag is considered equal to mx {n = 0,1,2,...) there 
will be no ratchet current, v, since in this case both the directions of x are identical in all 
respects. However, when (p y^ nn appreciable v is obtained in a small range of Fq with a 
peak in an intermediate Fq for given 70, T, and (p. 



In an earlier work 







the variation of ratchet current v as a function of the amplitude 



Fq of the applied square-wave forcing F{t), with a frequency of 5 x 10 cycles per unit 



time, was shown (Fig. 5 of Ref. 8|]) for the two cases, adiabatic drive and square- wave drive 
for 7o = 0.035 and temperature T = 0.4. The range of Fq over which ratchet current 
is obtained in the square-wave drive case was wider [0.1 < Fq < 0.8] compared to the 
adiabatic drive condition [0.07 < Fq < 0.12] and the peak current also occurs at a larger Fq 
value. The range, though wider, still remains well below F^, the critical field at which the 
potential barrier to motion just disappears. The current is, therefore, essentially aided by 
thermal noise. 

Since in the steady state ((Ax(t))^) ~ t we can define the diffusion constant D: 
((Aa;(t))^) = 2Dt. Interestingly, for given 70, T, and (f), the diffusion constant D shows 
nonmonotonic behaviour with the field amplitude Fq and it peaks around a value of Fq 
where v attains maximum at the drive frequency of 5 x 10~^ cycles per unit time. That is 
to say the ratchet current is maximised when the system is most diffusive. To compare the 
extent of this diffusive spread with the directional average displacement a quantity, Peclet 
number Pe, is defined as the ratio of square of mean displacement x in time t to half the 
square of diffusive spread in the same time interval t p^] : 

Pe= — = — . (1.1) 

Dt D ^ ' 

Thus, if Pg > 2 the motion is dominated by directional transport and hence it is considered 
coherent otherwise the net displacement is overwhelmed by diffusive motion. Our calculation 
shows that in the region where the ratchet current is appreciable and in particular where v 
peaks Pe is much larger than 2. Here the ratchet current is obtained when there was neither a 
bias to help the system nor any load to oppose it. The current, however, is not large enough 
and when a small load is applied against current the current either reduces to a small level 
or starts fiowing in the direction of the applied load. Thus, in the given circumstances, 
no appreciable useful work can be extracted from this inhomogeneous (frictional) ratchet. 
However, even in the absence of any external load the particle keeps moving against the 
frictional resistance. Leaving out the symmetric diffusive part of the motion the particle's 
unidirectional (ratchet) current v is maintained against the average frictional force. The ratio 
of this work (the ratchet performs against the frictional drag) to the total energy pumped 
into the system from the source of the external forcings is termed as Stokes efficiency, r]s, 
of the ratchet. 



An expression for r]s has been derived earlier |2, llO|, lUl which involves v as well as the 



second moment of the velocity v calculated from the probability distribution P{v) of the 
velocity v{t) recorded all through the trajectory of the particle. We have calculated rjs as 
a function of the amplitude of the applied forcing, rjs shows a peak, the position of which, 
however, does not coincide with the peak of the ratchet current. The distribution P{v) is 
almost symmetric about v = and the velocity dispersion grows monotonically approaching 
to be linear in Fq at large Fq. 

Recently, it has been reported [12] that in a tilted periodic potential an underdamped 
particle motion shows dispersionless behaviour in the intermediate time regime for a range 
of tilt values. In the present work we show that when the system is driven by a square- 
wave forcing of appropriate amplitude such dispersionless transient behaviour with added 
richness can be observed. The dispersionless behaviour of constant tilt gets punctuated 
and oscillatory behaviour of dispersion of different kinds, depending on the frequency of 
the periodic drive, naturally emerges. Interestingly, however, contrary to expectations, 
dispersionless particle motion do not contribute to (instead hinders) ratchet current in this 
system. 

In section II the basic equation of motion used in this model calculation will be presented. 
The section III will be devoted to the presentation of the detailed results of our numerical 
calculation. In the last section (Sec. IV) we shall conclude with a discussion. 

II. THE MODEL 

This part (II) of the work is an extension of our earlier work (part I), where the 
system was driven adiabatically 9|] to obtain ratchet current. In the present case we drive 
the system periodically by a symmetric square-wave forcing and calculate the particle 
trajectories x{t). The choice of square-wave forcing, instead of a sinusoidal forcing, is 
to make a direct contact with the adiabatically driven case. We thus have the forcing F{t) as, 

F{t) = ±Fo, {nT<t<{n+l)T), 

= tFo, ((n + i)r<t<(n + l)r), 

where r is the period of forcing and n = 0, 1,2, .... In what follows we shall refer to tra- 
jectories x{t) computed using the upper signs as the odd numbered trajectories and those 
calculated with the lower signs as the even numbered trajectories. 



The motion of a particle of mass m moving in a periodic potential V{x) = —Vosin{kx) in 
a medium with friction coefficient 7(x) = ^o{l — Xsin{kx + (/))) with < A < 1 and subjected 
to a square-wave forcing F{t) is described by the Langevin equation J13l|. 

Here T is the temperature in units of the Boltzmann constant ks- The Gaussian distributed 
fluctuating forces ^(t) satisfy the statistics: < ^{t) >= 0, and < ^it)^{t ) >= 26{t — t). For 
convenience, we write down Eq.(2.1) in dimensionless units by setting m = l, Vo = l, A; = l 
so that T = 2 corresponds to an energy equivalent equal to the potential barrier height at 
Fq = 0. The reduced Langevin equation, with reduced variables denoted again by the same 
symbols, is written now as 

S = "^^"^^^ + """'"^ + ^^^^ + ViWm, (2.2) 

where j{x) = 7o(l — Xsin{x + (j))). Thus the periodicity of the potential V{x) and also the 
friction coefficient 7 is 27r ij]. The potential barrier between any two consecutive wells of 
V{x) persists for all Fq < 1 and it just disappears at the critical field value Fq = Fc = 1. 
The noise variable, in the same symbol ^, satisfies exactly similar statistics as earlier. 

The Eq. (2.2) is solved numerically (with given initial conditions) to obtain the trajectory 
x{t) of the particle for various values of the parameters Fq, 70, and T. Also, the steady state 
mean velocity v of the particle is obtained as 

z;=(lim^), (2.3) 

where the average (...) is evaluated over many trajectories. The mean velocity is also calcu- 
lated from the distribution P{v) of velocities giving almost identical result. 

III. NUMERICAL RESULTS 

The Langevin equation (2.2) is solved numerically using two methods: 4*''-order Runge- 



Kutta[15| and Heun's method (for solving ordinary differential equations). We take a time 
step interval of 0.001 during which the fiuctuating force ^(t), obtained from a Gaussian 
distributed random number appropriate to the temperature T, is considered as constant 
and the equation solved as an initial value problem. In the next interval another random 



number is called to use as the value of ^ and the process repeated. By a careful observation 
of the individual trajectories of the particle shows that by t = 10^ the particle completely 
loses its memory of the initial condition it had started with. When we look for steady 
state solutions the trajectory is generally allowed to run for a time t ~ 10^. Therefore, for 
steady state evaluation for v, etc. the results become independent of initial conditions. The 
(Runge-Kutta) method had earlier been used and obtained correct results[l6| in a similar 
situation. Also, the method was checked against results obtained earlier for the adiabatic 
case (using „.tnx co„t,„ued f.act.o. method) a.d found to co,„pa.e we.. ^uaUtafvelyfl. 
Heun's method when applied in similar situations take much less time than the Runge-Kutta 
method and yields qualitatively as good result (Fig.l). With this confidence in our numerical 
procedures, we apply either one or the other of these two numerical schemes as the situation 
demands. We take A = 0.9, 70 = 0.035, and T = 0.4 all through our calculation in the 
following. 

The motion of the particle is governed by the applied square- wave forcing F{t). As F{t) 
changes periodically so does the position of the particle. In view of this effect we start our 
simulation at t = with F{0) = +|-fo|, and —\Fo\ for alternate trajectories (to be referred 
to respectively as odd numbered and even numbered ones, as mentioned earlier). This gives 
a nice nonoscillating variation of the overall average position when averaged over an large 
even number of trajectories. However, while calculating the position dispersions or velocity 
dispersions at a given time t, the even and odd trajectories are treated separately to calculate 
the deviations from their respective mean values. 

A. The ratchet current 

In Fig.l, v{Fo) are plotted for two values of phase lag (p = 0.3, and 0.4 for the same 
values of 70 and T to illustrate the effect of for r = 2000. The ratchet current v also 
shows nonmonotonic behaviour as a function of the period of the drive. In the inset (Fig. 
1) we plot the varation of t; as a function of the time period of the drive for Fq = 0.26. 
For these parameter values the current v peaks at period r ^ 1000. For comparison of 
time scales, it may be noted that for an equivalent RCSJ model of Josephson junctions the 
characteristic Josephson plasma frequency uj turns out to be about 10^ times larger than 
the drive frequency corresponding to T^ = 1000, where Tq = |. (The sign of Fq is changed 




FIG. 1: Shows the variation of v with Fq for phase lag (p = 0.3 and 0.4 with Tq = 1000. The inset 
shows the variation of v with time period Tq for Fq = 0.26. 

after every time interval Tq.) In this sense we obtain appreciable ratchet current only for 
very slow drives. It should, however, be noted that in the infinitely slow adiabatic case the 
ratchet current is effectively zero for Fq > 0.12 for 70 = 0.035 at T = 0.4. In the following we 
shall present the results obtained using (p = 0.35 except when mentioned otherwise explicitly. 



B. The steady-state dispersions 

The position dispersions ((Ax(t))^), where Ax{t) = x{t) — {x{t)) are evaluated over a 
large number of trajectories for various values of Fq, and Tq = 1000. It is found that the 
dispersions fit nicely to 

\og[{{Ax{t))')] = \og{t) +\og{2D), (3.1) 

for large t, typically t > 10^ (Fig. 6 of J8|). 

From the linear fit of the graphs we calculate the diffusion constants D{Fo) and the result 
is shown in Fig. 2. The diffusion constant has a large value between F ^ 0.15 and 0.35. The 
peak height is quite large ~ 800. As Fq is increased D decreases sharply and becomes smaller 
than 50 (which is less than 10% of its peak value) for Fq > 0.7. This [0.15 < Fq < 0.35] 
is also the region where the ratchet current v is appreciable. The Peclet number , P^, as 
defined earlier, are also calculated as a function of Fq. They are plotted in the inset of 
Fig. 2. It is clear from the figure that in the same region, Pg is also much larger than 2. This 
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FIG. 2: The variation of the diffussion constant D as a function of the driving amphtude Fq with 
Tf7 = 1000. The inset shows the variation of the corresponding Peclet number Pg with Fq. 



indicates that in the region [0.15 < Fq < 0.35] the particle motion is highly diffusive but 
concomitantly it is greatly coherent too. This is also indicated by the observation that even 
though the position dispersions (fluctuations) are large the relative fluctuations of position 
in this region are considerably low (< 1). As indicated by the result in the adiabatic case 
(Fig. 3 [8]) this range of Fq of coherent motion is expected to shift as the value of 70 is 
changed. 

Though our system is different from that of Machura, et. al. [2|, at this point it would 
be interesting to make a comparison with their result. They observe that for their low 
temperature case Dq = 0.01 in the vicinity of a ^ 0.6 the velocity fluctuation underwent 
a rapid change (Fig. la of [2]). To translate this to our case[lj] a ~ 0.6 is equivalent to 
Fq ~ 0.2 and given their potential barrier being just about half of the value in our case 
one should expect the peaking of velocity dispersion to occur below Fq = 0.4. Taking into 
consideration of our temperature (T = 0.4) being 40 times 0.01 the phenomena should occur 
much below Fq = 0.4. In this sense the region [0.15 < Fq < 0.35] seems quite reasonable. 
Also, V of Fig.3a of Ref. 2l] at Dq = 0.4 make a good comparison with Fig.l in our case. 
However, as mentioned earlier the two systems are quite different in basics to have an exact 
comparison. 

The velocity distribution P{v) also shows interesting behaviour. In Fig. 3, we plot P{v) for 
three values of Fq. A sharp peak which is almost indistinguishable from a Gaussian centred 




FIG. 3: Plot of velocity distribution P{v) for three values of driving amplitudes Fq = 0.05, 0.12 
and 0.30. The figure in inset shows the variance of velocities as a function of Fq fitted with a 
straight line to show the linear growth of variance at large Fq . 

at v = for small Fq = 0.05 gets split up into three peaks for Fq = 0.12, and similarly for 
Fq = 0.30, with the central peak, gradually diminishing. This shows a behaviour, including 
the nearly linear growth of the variance with Fq (inset. Fig. 3), quite similar to what has been 
reported earlier in a different system|2]. There is, however, one difference. The side peaks of 
P{v) in our calculation have origin in the running states of the particle. It is, perhaps, due 
to the square-wave drive, instead of sinusoidal drive, that for as low amplitude as Fq = 0.3 
we get three disjoint velocity bands and at Fq = 0.6 we get just two bands, the central band 
being almost unpopulated. The three peaks, for example for Fq = 0.3, could be fitted to a 
combination of three Guassians. With a cursory look, the left and right Gaussians barely 
show much difference. However, 



(v) 



vP(v)dv 



(3.2) 



gives approximately the same value as v, and {v){Fq) showing exactly the same nature as 
v{Fq) (Fig.4). 
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FIG. 4: Shows the variation of the steady state mean velocity f , Eq. (2.3) and (v), Eq. (3.2) for 
the same parameter values as in Fig. 2. 




FIG. 5: Shows Stokes efficiency, rjs as a function of Fq for the same paramater values as in 
Fig. 2. The inset shows the difference in the velocity distribution for symmetric (three peaks) and 
asymmetric drive for the same value of Fq = 0.16 and r = 2000 and a = 0.2. 



C. The efficiency of ratchet performance 

From the velocity distribution P{v) we calculate the Stokes efficiency, tjs, defined as 

/.,\2 



Vs 



\{v')-T\' 



a, 

(3.3) 
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as a function of Fq. Fig. 5 shows that rjs is larger in the same range of Fq where it shows 
larger v. The peak of r]s, however, does not occur at the same position as the peak of v. 
It is, however, to be noted that the plotted figure is calculated from averages over a small 
number (~ 20) of ensembles because it is computationally quite expensive to obtain results 
for the steady state (maximum t = 10^) and hence not feasible to obtain averaging over a 
larger number of ensembles. Though the qualitative behaviour is encouraging the efficiencies 
are small ~ 10~^. In the adiabatic drive case (part I, Ref. 9|) we have found that Stokes 
efficiency depends on various parameter values: 70, T, etc. The efficiency shown here is 
for a small 70 = 0.035, Tq = 1000, and T = 0.4 where the current is also very low. The 
efficiency of this symmetrically driven system can, however, be improved to a good extent 
by an optimal choice of these parameters. 

An inertial ratchet driven by a zero mean asymmetric drive can, however, give a highly 
efficient performance compared to the symmetrically driven ratchet. For example, when the 
system is driven by a field 

F{t) = ±Fo, {nr <t<{n + a)r), 

= T^Fo, {{n + a)T<t<{n+ l)r), 

with a = 0.2 gives an efficiency of 3.8 x 10^^ compared to 6.2 x 10^^ in the symmetric-drive 
{a = 0) case with Fq = 0.16 and r = 2000. This is made possible because in the symmetric 
drive case the particles move on either direction with almost equal probability whereas in 
the asymmetric drive case the particle motion in one direction is practically blocked, as 
is evident from the corresponding velocity distributions shown in the inset of Fig.5. The 
contribution of the system inhomogeneity for this improved performance is, however, quite 
insignificant. 

D. The transient-state dispersions and the ratchet current 

When a constant force F is applied to the system it shows dispersionless behaviour: 
{{Ax{t))'^) does not change with time in the intermediate time scales, roughly [10^ < t < 10^], 
for around [0.12 < F < 0.7] at T = 0.4 for 70 = 0.035. The result of dispersionless behaviour 
had originally been shown and explained |12| beautifully for constant friction 7 case: the 



position distribution moves undistorted at constant velocity v = — or equivalently, velocity 
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FIG. 6: The plot of position dispersions ((Aa:(t))^) versus time t (in logarithmic scale) for different 
values of (Tq) of forcing with Fq = 0.2. The inset shows the clipped part of the plot at larger time. 



distribution remains undistorted centerd at f = — . The interval [ti < t < t2] of time t 
during which the system shows this remarkable intermediate-time behaviour depends on 
the tilt force F, as should also on other parameters, ti is roughly of the order of but much 
larger than the Kramers passage time corresponding to the lower of the potential barriers on 
either side of a well. The transient-time dispersionless particle-motion behaviour is sensitive 
to initial conditions. In the following we specifically begin from the bottom of the well at 
x = I with particle velocities appropriate to the Boltzmann distribution at temperature 
T = 0.4. 

When the inhomogeneous system is driven periodically by a sqaure-wave forcing of am- 
plitude Fo, the dispersionless coherent nature of average motion gets interrupted depending 
on the period Tq of the forcing [Fig.6]. When ti < Tj^ < t2, at t = Tq the dispersion gets a 
jerk and shoots up only to get flattened again to an another bout of dispersionless regime. 
This regime too gets a similar jolt after another Tq and the process continues for a large 
number of periods. When the direction of the applied force is changed the 'forward moving' 
particles are forced to halt momentarily to begin moving in the new direction of the force 
afresh. While in the state of halt particles are more likely to find themselves closer to the 
bottom of some well and thus the system gets initialised as in the beginning. The system 
finds itself in similar situation again and again periodically with each change of force direc- 
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tion and continues with its unfinished dispersionless sojourn for a large number of periods 
with remarkable robustness [Inset of Fig.6]. However, when Tq < ti the system never gets 
a chance to experience its dispersionless journey because only a fraction of the particles 
get the opportunity to acquire the required constant average velocity[12.] of — and the rest 
keep lagging behind even by the end of constant force duration Tq. Instead, as soon as the 
direction of the force F is reversed, after Tq the dispersion dips after a brief climb up, as 
the particles get herded together briefly before getting dispersed further in the reversed di- 
rection of Fq. This can be seen very clearly in the time evolution of the position probability 
distribution profile P{x,t). The front of the P(x,t) moves with velocity — while the rest 
lag behind it moving at a slower speed but trying to catch up with the front throughout 
Tq. This process of dispersion dipping (after a small contunuing rise) and rising to a higher 
value after each Tq is repeated for several tens of periods [Fig.6]. 

The intermediate-time dispersionless motion is not an exclusive characteristic feature 
of inhomogeneous systems. It is a characteristic feature of inertial washboard potential 



system 12|]. However, its study in the inertial inhomogeneous system provides a convincing 
explanation of the variation of ratchet current as a function of Tq [inset of Fig.l] and helps 
in finding a criterion to improve the performance of the ratchet. 

In fig. 7 the average displacement of particles as a function of time when driven by equal 
number of ±i^(t) profiles is presented for Fq = 0.2, and T^ = 5000. This case corresponds 
to the repeated dispersionless motion shown in Fig.6. Fig. 7 clearly shows that during the 
dispersionless motion the average dispacement of particles effectively remains constant. In 
other words, during the period of dispersionless motion the particles move equally in the 
left as well as in the right direction thereby contributing nothing to the ratchet current: 
while in the dispersionless motion the particles fail to see the frictional inhomogeneity of 
the system. All the change in the average displacement and hence all the contribution to 
the ratchet current comes during the dispersive period of motion. This is shown in the inset 
of Fig. 7 where for clarity the mean particle positions for F{t) beginning with F{0) = —Fq 
(even numbered trajectories) are shown as a function of time with their sign reversed. The 
mean particle displacements for odd and even numbered trajectories differ only during the 
interval just after the reversal of Fq and before the dispersionless regime begins and the two 
lines of mean positions (insets of Fig. 7) run parallel during the dispersionless regime. This 
clearly indicates that in order to get a larger current an optimum choice of Tq needs to be 
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FIG. 7: The average displacement of particles as a function of time, driven by equal number 
of iF{t) profiles (or equal number of odd and even numbered trajectories) for Fq = 0.2, and 
Tq = 5000. The insets highlight the contributions to the mean displacement of odd and even 
numbered trajectories separately, leading to the main figure. The mean displacements for the even 
numbered trajectories are shown with a reversed sign. 



made which, naturally, avoids the dispersionless regime but is not too small in order to allow 
the particles to leave their potential wells. This conclusion is well supported by the inset of 
Fig.l. 

The velocity dispersions and position dispersions together show interesting behaviour. 
Fig. 8 shows that during the dispersionless regime when the position dispersion is constant 
and maximum the velocity dispersion is also constant but it has a minimum value. This 
minimum constant value is repeated in all the hTq, n = 1, 2, ... intervals whereas the value of 
the constant position dispersion increases in every successive uTq interval as shown in Fig.6. 
In the dispersive regimes the velocity dispersions are squeezed to very sharp troughs exactly 
where the position dispersions show sharp peaking. In the inset of Fig. 8 these dispersions 
are shown for Tq = 250. The onward rush of the particles do not halt immediately after 
the direction of Fq is changed at uTq but it continues for a very short time giving a small 
increase in the spread of P{x). Then a majority of particles stop, giving a sharp peak in the 
P{v) at i; = reducing its spread drastically. At that moment the product of position and 
velocity distribution spread becomes a minimum. The reverse journey thereafter increases 
the spread of P{v) but there is a slow squeezing of P{x) before it begins to spread again. The 
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FIG. 8: Illustration of velocity dispersions ((At;(t))^) (thin line) and position dispersions ((Ax(t))^) 
(bold line) during a time interval for Tq = 5000. The inset shows the corresponding plots for square 
drive forcing with smaller Tq=250 with no dispersionless regime. The thin vertical lines roughly 
indicate the positions of the extrema of the dispersion curves. 



maximum P{x) squeezing, however, does not exactly coincide with the largest of the broad 
P{v) but it is at a rather closer range. In this case too the minimum velocity dispersion 
remains constant for all uTq. But the wings of P{x), though thin, keep spreading with time 
giving an average increase of dispersion as time increases. However, most of the particles 



l^ol' 



1^0 1 



remain confined roughly to a region [—^—^Tq < x < +—^Tq] for a long time. 



IV. DISCUSSION AND CONCLUSION 

The ratchet effect, in this work, is brought about just by the phase lag (p between the 
potential and the friction of the medium, without having to have an external bias. This is 
seemingly a weak cause to generate unidirectional current. The Figs.l through 5 refer to a 
square- wave forcing with T^ = 1000. The choice of this Tq clearly avoids the dispersionless 
regime. Yet, this is not the optimum value of Tq. It should have been around 500 in order 
to get the largest possible ratchet current. This choice would have definitely enhanced the 
efficiency of operation. The same can also be said about other parameters, such as T, and 
(p for 7o = 0.035. However, with the help of these figures we have been able to exhibit the 
qualitative trends shown by the ratchet. 
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In the inset of Fig. 6 we have drawn a straight hne with slope 1 as a guide to show that 
ultimately the curves should achieve that average slope at large times for various Tq values of 
drives. Even though the average slope of the curves have not yet reached the diffusive slope 
of one the small Tq^ curves are slowly approaching that value. One can,therefore, safely infer 
that in the steady state situation the effective diffusion constant should increase monotoni- 
cally with Tq for small Tq. The frequency of drive or equivalently Tq, thus, plays important 
role about how the particles diffuse out of their wells. For example, The population of the 
initial well depletes with time exponentially, N{t) = iV(0)e"''*, with b = 0.0023 for Tq = 250 
and b = 0.002 for Tq = 500, for 70 = 0.035 at T = 0.4 that we have studied. By the time the 
well gets effectively exhausted the first particles would have moved farther than a thousand 
of potential wells. Of course, this first well itself (as all others) keeps getting repopulated 
all the time. 

The dispersive behaviour for drives with Tq > t2 is difficult to study because it takes a 
very large computer time to arrive at a concrete result. However, the indications are there 
that for these large Tq also, the system will show repeated dispersionless regimes, though 
somewhat enfeebled because the process of diffusion will dominate at these large times. 

To conclude, the study suggests an interesting method of obtaining ratchet current in 
inertial noisy systems by exploiting the frictional inhomogeneity of the medium. It also ex- 
hibits clearly that the ratchet current, in this system, is contributed by dispersive conditions 
and not by coherent movements of particles. 
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